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Abstract. We investigate the interplay of quark and meson degrees of freedom in a physical state repre- 
senting a near-threshold resonance for the case of a single continuum channel. We demonstrate that such a 
near-threshold resonance may possess quite peculiar properties if both quark and meson dynamics generate 
weakly coupled near-threshold poles in the S-matrix. In particular, the scattering t-matrix may possess 
zeros in this case. We also discuss possible implications for production reactions as well as studies within 
lattice QCD. 

PACS. 12.38.Lg - 11.55.Bq - 12.39.Mk 



1 Introduction 

With recent developments of -B-factories, a wealth of new 
charmonia states was reported, with properties of some of 
those being incompatible with simple quark model predic- 
tions. Exotic explanations, such as hybrids or tetraquarks, 
are suggested to describe the new states. On the other 
hand, these new enigmatic states are above open charm 
threshold, so that the spectrum of the "quenched" quark 
model is to be significantly distorted by charmed-meson 
pairs (for reviews, see e.g. Refs. [Tl|2j|3l|4l|5] ) • 

Moreover, some of the new states reside in the vicin- 
ity of 5*- wave thresholds. A most prominent example here 
is the famous X{3872) state, first observed in i3-meson 
decays [B]. The X(_3872) is extremely close to the DD* 
threshold, with DD* being in the S'-wave, if the quan- 
tum numbers of the X are l"*"^, as suggested by the data. 
There are also the y(4260) and y(4325) vector states 
[IlIH], with relevant_ S'-wave thresholds being the DDi at 
4.285 GeV and D*Do at about 4.360 GeV (Di and Dq are 
correspondingly 1+ and 0"*" D- mesons). The purely exotic 
charged state Z+(4430) [9 is close to the D*Di thresh- 
old and, for the charged Z+(4050) and_Z2+(4250) ^^the 
relevant thresholds are D*D* and DiD, respectiveljo- 

Threshold proximity implies that, independently of the 
binding mechanism, there should be a significant compo- 
nent of a hadronic molecule in the wave function of the 
state. So the question arises of how to distinguish between 
a genuine "elementary" particle (gg, hybrid, or compact 



^ These charged states are to be considered with caution, 
as seen only in one experiment (Belle). Besides, as quantum 
numbers of the Z-particles are not known, the thresholds are 
not necessarily the S-wave ones. 



tetraquark) and a composite state (hadronic molecule), 
and how to estimate the admixture of the latter. It was 
suggested in Ref. [TTlfT^ITB] that it is possible, in the case 
of a near-threshold bound state, to answer this question in 
a model-independent way: the state is mostly elementary 
if the effective radius is large and negative. The approach 
was generalised in Ref. P^[T3] to the case of presence 
of inelastic channels, as well as to the case of an above- 
threshold resonance. Related to this is the pole counting 
approach ^16] , in which the structure of the near-threshold 
singularities of the scattering amplitude is studied. It ap- 
pears that the state is mostly elementary if there are two 
nearby poles in the scattering amplitude, while a compos- 
ite particle corresponds to a single near-threshold pole. 

The approaches of Ref. [TTl[T2l[T3l[T4l[T5l[T6] are based 
on the effective-range expansion of the scattering ampli- 
tude and, as such, are expected to be valid for the mo- 
menta involved much smaller than the inverse range of the 
force. However, as noticed in Ref. [TTllT^fT^ . the effective- 
range formulae are not the most general ones even in the 
small-momenta limit. The scattering amplitude, as a func- 
tion of energy, can have a zero and, if this zero is situated 
in the near-threshold region, the effective-range expansion 
would fail. In the present paper we develop a formalism 
which allows one to pinpoint the source of this failure, to 
study the hadronic observables in the presence of this zero, 
and to identify the physical situation in which it occurs. 



2 General formalism 

We consider a physical state which is a mixture of a bare 
state {qq or compact tetraquark) and a dynamical (mole- 
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cule) component and represent its wave function as: 



cIV'o) 
x\MiM2 



(1) 



where lipo) is the bare elementary state with the probabil- 
ity amplitude c, while xip) describes the relative motion 
in the system of two mesons (M1M2), with the masses mi 
and 7712, respectively, and with the relative momentum p. 
The wave function \^) obeys a Schrodinger-like equation: 



H\^)=E\^), 



with the Hamiltonian 



(Ho Vqh 



where 



mi + m2 + E{) is the bare state mass, and 
2 

2/i 



(2) 

(3) 
(4) 

(5) 



where /x is the reduced mass. 

The term Vqh is responsible for the dressing of the bare 
state, which is given by the transition form factor f{p): 



(V'olVlMiMa) = /(p). 



(6) 



The Schrodinger-like equation ^ is equivalent to the 
system of equations for the c and x, which reads: 



c{E)Eo + / I{p)x{p)d'p - c{E)E, 
^Xip) + <E)}{p) + j V{p, k)x{k)d\ = Exip). 

^ . . . 

On substituting c{E) from the first equation to the 

second one, we arrive at the Schrodinger equation in the 
mesonic channel, with an effective potential 

VMp,p',E) = V{p,p') + 4?^. (8) 

The off-shell mesonic ^-matrix t{p^',E) is a solution of 
the Lippmann-Schwinger equatiorp. 



t{p,p',E)^VMp,p',E)~ J 



V,s{P,q,E)t{q,p\E) 
q'^/{2^i) - E -iO 

(9) 



The solution of this equation can be written as 

+ (10) 



Normalisation of the t-matrix is such that the M1M2 scat- 
tering amplitude is given by f{k,k,E) = -4tt'^ fj,t{k, k, E), 
with E = P/(2/i). 



where tv{p,p' , E) is the t-niatrix for the potential prob- 
lem. 



tvip,p',E)^Vip.,p')~ J d\ 



V{p,q)tv{q,p',E) 
q^/{2^) - E -iO 

while the dressed vertex functions are 



{p,E)^f{p)~ j 



tv{p,q,E)f{q) 
-q'^/i2fi)-E-iO' 



and 



g{E) = [ d\- 



q^/{2n)-E-iO' 



, (11) 



(12) 



(13) 



(14) 



'q^/{2fi)-E-iO' 

The system of equations (O can possess bound states 
(generally, more than one). For a bound state i with the 

binding energy e^-* , the solution of the system ([7]) may be 
written as 



cosd,, Xs (P) = ^i(p)sin6'-o (15) 



where ^i{p) is normalised to unity. Then the wave func- 
tion of this bound state, 



cos 9i\i^o) 
sine, ^Pdp)\MiM2)J ' 



is also normalised, = 1- 



The quantity 



(17) 



introduced in Ref. pTlfT^fTSi , gives the probability to find 
a bare state in the wave function of the bound state i. 
The solution of the system ([7]) for the continuum is 



Ck{E) = -^-1^ J d^pxk{p).f{p) = 



{k,E) 



E-Ea + g{Ey 



t{p,k,E) 



(18) 



fc2 

p2/(2^)-S-iO' ^"2jl' ^^^^ 



Xk{p) = 5{p-k) - 



so that the continuum counterpart of the quantity (jl7p . 
the spectral density w{E), gives the probability to find the 
bare state in the continuum wave function |17| . It can be 
found as: 

w{E) ^4TTfik\ck{E)\^e{E), (20) 

with Ck{E) given by Eq. ^E^. As shown in Ref. [T7], the 
normalisation condition for the distribution w{E) reads: 



/"OO 

/ w{E)dE ^l-J^Z^^ 
Jo 



(21) 



where the sum goes over all bound states present in the 
system, while the corresponding Z-factors are given by 
Eq. dni). 

We are interested in the energy range close to the 
threshold. To perform a low-energy reduction, let us as- 
sume that the scattering length approximation is valid for 
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the potential problem, so that the potential f-matrix ty 
takes the form: 



1 



47r2/i(— fly — ik) 



(22) 



where ay is the scattering length for the potential V{p,p') 
and ellipsis stands for the range of forces corrections. We 
use the sign convention, where a negative scattering length 
corresponds to an attractive potential, however, too weak 
to produce a bound state while a positive scattering length 
reverses to either a repulsive potential (which cannot pro- 
duce a singularity of the S'-matrix, thus we automatically 
have 1/ay ~ /3, /3 is the range of forces) or a bound state, 
which is the case of interest here. Note, the above expres- 
sion is useful only for ay/3 3> 1 and then applicable only 
for momenta of order 1/ay. In this case ty has a near- 
threshold pole. However, it should be stressed that this 
singularity is observable only in the weak coupling limit 
to the quark states located nearby, for otherwise all singu- 
larities undergo a significant mixing, the origin of which 
is the appearance of iy in the dressed vertex functions of 
Eqs. and (fT5|) . For a detailed discussion of this effect 
see Ref. [18l[l9]. 

It is convenient to define the loop functions: 



g{E) = j d\ 
g'[E) = j d\ 



I [2^1) - E ^ iQ 

fjq) 

gV(2/x) -E-iO 



/2(i? + 47rVfc), 



(23) 



= /o(i?' + 47rV«fc), 



where the constants R and R' are expected to be of order 
/Lt/3 and corrections 0{k'^/f3'^) where dropped; in addition 
we introduced /o = /(O). Defining 



Rv 



ay ' 



(24) 



one arrives at the following form for the on-shell ^-matrix 



E^Ec 



[E - Ea]Rv + fi[RRv - R''^] + '^^:^^^lik[E - Ec] ' 

(25) 

Ec — Eq — fQ{R + Ry — 2R ). 

Expression (|25p for the scattering amplitude has a zero at 
E = Ec, and the effective-range expansion for the am- 
plitude is not valid, if Ec is in the near-threshold 
region. 



3 Flatte formulae 

In this chapter we develop a generalisation of the well- 
known Flatte parameterisation [5n| for the near-threshold 
scattering amplitude, and express the i-matrix of Eq. (j25p 
in terms of Flatte-type parameters. 

Let us define the Flatte energy Ef such that the real 
part of the denominator in Eq. (j25D has a zero when E = 

[Ef-Eo]Rv + fo[RRv-R'^] = 0, (26) 



and get rid of Eq this way. Then one can express Ec in 
terms Ef. 



f2 

Ec = Ef - -5— (^' - Rv) , 



and, using the identity 

E — Ef E — Ef 



{E-Eff 



E — Ec Ef — Ec {E — Ec){Ef ~ Ec) 
rewrite the f-matrix in the Flatte-type form as: 



t{E) = 



9f 



Sn^fiVpiE)' 



Here 



Ef 



Ec 47r2^/2 



Ry 



and 



Vf{E) ^E -E 



f 



Rl 



(E-E 



{R — R\, 



f) 



E-Er 



(27) 



(28) 



(29) 



(30) 



(31) 



If \Ec\ ^ \Ef\, we are back to the standard Flatte ap- 
proximation for the near-threshold scattering amplitude 



tF{E) 



1 



9f 



Sn'^liE-Ef + ^gfk' 



(32) 



which depends on two parameters, gf and Ef, in contrast 
to the expression ([^. which depends on three parame- 
ters: gf, Ef, and Ec- One can see that, for \Ec\ ^ \Ef\, 
the standard Flatte formula is severely distorted by the 
presence of the i-niatrix zero, with a drastic effect on the 
behaviour of the elastic scattering cross section. 

The Flatte-type formulae given above define the t- 
matrix in terms of the parameters Ef, gf, and Ec- Fol- 
lowing S. Weinberg pT[|T?llTB] . it is instructive to consider 
the case of a single bound state present, and to express the 
t- matrix in terms of the binding energy es, the Z- factor 
(which is the probability to find a bare state in the physi- 
cal bound state wave function), and the scattering length 
ay. 

We start from the expression for the inverse i-matrix 
(EH), 

E ~ E 

t-^{E) = -J-Rv + 47rVifc- (33) 



E-Ec 

In the vicinity of the bound-state pole one has 



t-\E) 

with (see Ref. [HIIIIIIS]) 



es + E 

7 1 



This defines 



Z 



2 VS^es,, „^ 
2eB 



1-Z EB+Ec 



1 - 



iv 



^2^eB J ' 



(34) 



(35) 



(36) 
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with 7v = l/fly being the inverse scattering length for 
the potential problem in the absence of the quark state. 
Thus the following formulae can be obtained [11,12,13 : 



k cot 5 

and 

Ec = 



y^2jIe^{E + eB){eB + E, 



■c) 



-es 1 



2eBiE~Ec) 
2(1 -Z) , 2(1 -Z) 7y 



Z 



1 - Z' 
(37) 



(38) 



Now, for \Ec'\ ^ eb, the effective-range expansion is 
recovered from Eq. ([57]) , with the scattering length a and 
the effective radius given by: 



2(1 -Z) 1 
{2-Z) ,/2i^' 



{l-Z) ^2^ 



(39) 



which allows one to define the quantity Z in a model- 
independent way, if the effective-range parameters a and 
Te are available. 

On the contrary, for \Ec\ < es the effective-range 
expansion does not converge anymore for energies larger 
than \Ec\- Correspondingly, formulae from Eq. do 
not hold anymore. In the next section we discuss the cir- 
cumstances, when the described unusual behaviour may 
occur. 



4 Pole structure in the presence of the 
t-matrix zero 



Due to a very simple relation (j30p between Ec and 7v, 
it is instructive to study the amplitude singularities in 
terms of the variables Ef, gf, and 7y. The expression for 
the scattering t-matrix in this case takes the form: 



t 



1 



47rV [E - Ef){-fv + ik) + ^gnvk' 



(40) 



As seen from expression PO)) . there generally could be up 
to three near-threshold poles, and the equation defining 
the pole positions in the fc-plane reads: 



(fc^ — 2fj.Ef){'-fv + ik) + ifigjjvk = 0. 



(41) 



Denote the solutions of Eq. (|^T|) as fci , ^2 , and ^3 • As 
follows from the explicit form of Eq. (|4T|) , one of these so- 
lutions is always imaginary, while the other two are either 
both imaginary or placed at the lower half-plane symmet- 
rically with respect to the imaginary axis, as required by 
the general properties of the S'-matrix. Thus, the set of 
parameters {Ef,gf,'-fv} used in Eq. (|iD)) is fully defined 
by the poles of the scattering i-matrix {ki, ^2, fcs}: 



9f = 
7y = 



1 



kik2k-i 



2fj, ki + k2 + ks 

(fcl +fc2)(fci +fc3)(fc2 +fc3) 

i^{ki + k2 + ^3)2 
-i{ki + k2 + fca). 



(42) 

(43) 
(44) 



The position Ec of the t-niatrix zero can be also ex- 
pressed in terms of {fci, ^2, fca}: 



Ec — Ef 



1 



-^{kik2 + kik3 + k2k3). (45) 
2/i 



Notice that, for the given Ef and -fv, the range of values 
for Ec is restricted by the condition g/ > (see Eq. ([50]) ). 
which follows from causality. 

It is illustrating to discuss the behaviour of the three 
poles when the effective coupling 5/ is varied from zero to 
some large value. To study the properties of the system 
we introduce the scale A <^B that defines the region of 
applicability of the equation^ If the coupling of the bare 
state to the mesonic channel is switched off, gf — 0, the 
solutions of Eq. (|^T|) are 



^3°^ = ijv- 



(46) 



The first pair of poles corresponds to the bare quark state 
(uncoupled from the mesonic channel) with the energy 
Ef = Eq, and the third one is due to the direct interac- 
tion V{p,p') in the mesonic channel, as already described 
above. Note, however, that the numerator of the i-matrix 

(see Eq. (00])) has a zero exactly at the momenta k^fl 
which eliminates the contribution of the corresponding 
poles. Thus, in the limit of an extremely small coupling, 

<C |i^/|/|7y| (see Eq. (05])), there is only one pole k^ in 
the scattering t-matrix that corresponds to a bound or vir- 
tual state in the mesonic channel, depending on the sign 
of 7y. As the interaction g/ is switched on and increases, 
gf \Ef\/\jv\, the poles start to move and to mix, and 
the t-matrix may have a zero and three poles in different 
places of the fc-plane. Note also that in this case there is 
no straightforward interpretation of 7y possible anymore 
due to the mixing of poles. Clearly, the zero of the t-matrix 
can have an impact on the observables only if it is in the 
region of applicability defined as A. From Eq. (|^5]) one can 
see then that Ec ~ A only if all three poles are located 
very close to the threshold. In other words, if either the 
potential pole or the quark poles are outside the region of 
applicability of our formalism, also the t-matrix zero lies 
outside that region. In the latter case the effective-range 
approximation is adequate for the problem under consid- 
eration, so that one can set Ec ^ E, A in all equations, 
thus recovering the standard Flatte formulae. 

Let us now assume that all three bare poles appear 
in the vicinity of the threshold {\Ef\ ~ 7y/(2/i) < A) 
but the coupling is strong {gf » |£'/|/|7y|). Clearly, in 
this case Ec ^ A and we again recover the standard 
lineshapes together with the standard effective range ex- 
pansion; only a single pole remains in the near-threshold 
regime — the expressions of Eq. (15^ hold again and we 
find, as expected, that this single pole corresponds to a 
state of predominantly dynamical (^molecular) origin. 

To summarise, the most peculiar situation takes place 

if: 



Note, however, that this applicability region might be fur- 
ther restricted by additional scales present in a particular prob- 
lem under consideration. 
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1. the bare quark and potential poles appear acciden- 
tally close to each other and to the threshold {\Ef\ ~ 
7^/(2m)<Z\); 

2. the quark state is relatively weakly coupled to the 
hadronic channel {gj ~ |£'/|/l7y|)- 

Thus, the appearance of a zero in the scattering amplitude 
very close to threshold, signaled also by an early break- 
down of the effective-range expansion, is a clear signal of 
the presence of both potential as well as quark poles. 

To detail these statements and to mimic t he pecu liar 
situation described above we choose |7y| ~ y^2fi\Ef \ for 
gf = 0. We choose for the reduced mass = 1 GeV as 
is appropriate for the X{3872) case, and define the near- 
threshold region by Z\ '-^ 1 MeV — the energy range rele- 
vant for the X{3872) particle (in case of the X(3872) the 
range of validity of the approach is not set by the range of 
forces, but by the closet threshold, which is only 8 MeV 
away). We then study the system for different values of 
gf for six representative scenarios (recall: 7y < refers to 
an attractive potential, which does not support a bound 
state, while 7y > refers to a bound state present from 
the pure potential scattering). For each case we fix the 
parameters Ef and jv, and change gj from very small 
values to the ones for which the t-matrix zero leaves the 
near-threshold region. 

- Case (i): Ej > 0, -fv > {Ef = 1 MeV, 7y 
45 MeV). See Fig.HKl 

- Case (ii): £;/ > 0, 7y < {Ef = 1 MeV, 7y = 
-45 MeV). See Fig. O 

- Case (iii-a): Ef < 0, 7y > 0, yj2^\Ef \ > {Ef = 
-1 MeV, 7v = 40 MeV). See Fig . Bl 

- Case (iii-b): Ef < 0, 7y > 0, ^j2[i\Ef \ < 7y {Ef = 
-1 MeV, 7y = 50 MeV). See Fi g. El 

- Case (iv-a): Ef < 0, 7y < 0, y/2n\Ef\ > \jv\ {Ef = 
-1 MeV, 7y = -20 MeV). See Fig. El 

- Case (iv-b): Ef < 0, jv < 0, ^/2fi\Ef\ < \jv\ {Ef = 
-1 MeV, -fv = -55 MeV). See Fig. El 

For cases (i) and (ii), the hadronic cross section van- 
ishes at E — Ec accompanied by a very peculiar energy 
dependence. For small couplings, the spectral density as 
well as the elastic cross section peaks at around Ef and 
the bound state, present in case (i), is a purely mesonic 
one. Then, with the increase of the coupling, the afore- 
mentioned bound state leaves the near-threshold region 
and acquires a large admixture of the quark component. 
In this strong-coupling regime the dynamics is defined by 
the presence of a single pole which corresponds to the vir- 
tual state. Due to the normalisation condition (I^TI) . the 
spectral density decreases with the increase of the cou- 
pling. 

It is instructive to discuss, for example. Case (i) in 
some more detail. In particular, in Fig. [2l the movement 
of the poles in the s-plane is depicted on both first (left 
plot) and second (right plot) Riemann sheets. The arrows 
indicate how the poles start to move when 5/ is increased. 
Clearly this plot is equivalent to the pole movement in the 
/c-plane depicted in the upper left panel of Fig.[Tl As usual, 
for a narrow resonance (small gf) only the pole in the lower 



half plane on the second sheet (see Fig. [5]) influences the 
physics. Here this state is nearly a pure quark state. As gf 
gets larger the width of the state grows (and the mesonic 
admixture increases). At latest, when the real part of the 
pole position reaches the threshold (indicated in the figure 
as the perpendicular, dashed line), both poles are equally 
important. Notice that, in this regime, the interpretation 
of the imaginary part of the pole position as half of the 
width of the state is lost. Similar pole movements in the 
subthreshold regime were reported previously in Refs. |21l 

In Cases (iii-a,b) one has Ec < 0, so that the i-matrix 
zero does not manifest itself in the hadronic cross sec- 
tion defined for positive energies only. In the meantime, 
with \Ec\ ~ ^, there are two bound states in the near- 
threshold region. With the increase of the coupling, only 
one bound state survives in the near-threshold region, 
with Z 0. The spectral density is small for all values of 
the coupling, as two bound states saturate the normalisa- 
tion condition (HI]). 

Finally, in cases (iv-a) and (iv-b) and for small cou- 
plings, the bound state is almost purely of quark nature, 
and Ec < 0. With the increase of the coupling, Ec be- 
comes positive and shows up in the hadronic cross sec- 
tion, while the bound state acquires a significant mesonic 
admixture. Further increase of the coupling forces the t- 
matrix zero to leave the near-threshold region, and the dy- 
namics is defined by a single bound state, which is purely 
molecular, with Z — > 0. 

So, for all cases considered, there is only one near- 
threshold pole in the strong-coupling regime which cor- 
responds to the scattering-length approximation for the 
mesonic t-matrix, independently of the underlying dynam- 
ics. 



5 Production reactions 

Unfortunately there is no experimental possibility to study 
elastic scattering of, say, charmed mesons, and our knowl- 
edge of the resonance properties comes from production 
experiments. In general there are two production mech- 
anisms possible, namely a production of the resonance 
of interest via its hadronic component or via its quark 
component. If the production source can be considered 
as point-like, then the production amplitude of the me- 
son pair (M1M2) through the resonance from the former 
mechanism can be written as 



Mh{E)^Th 1 



t{p,k,E) 



p^/{2fi) -E-iO 



|fcV(2M)=-B 

where J-h is the initial state production amplitude from 
the point-like source, unity stands for the Born term, and 
the second term defines the final-state interaction — see 
the left panel of Fig. [51 For a point-like source and small 
energies, the amplitude (|Tfl) is 



Mh{E)=Th{l-L{E)t{E)), 



(48) 
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L being the loop function describing propagation of the 
intermediate mesonic state, 



one finds: 



L{E) = 47rV('o + jfc), 



(49) 



with Iq ~ 13, where f3 is the range of the force. We may 
therefore write 



{E-Ef){Ef-Ec)-\logf{E-Ec) 
{E - Ej){Ef - Ec) + ^kgf{E - Ec) 



(50) 

Thus, in Mh{E) the zero in the production amplitude is 
shifted with respect to the zero of the f-matrix. However, 
for the energies Ef ~ Ec ^ /? this shift is small. The sec- 
ond term in the numerator of Eq. (1501) dominates, and the 
production rate through the hadronic component alone is 

— = const X k\t{E)\ 0{E) 

dE 

where the denominator Dp is given by Eq. (j3ip . Thus 
one might expect that all said above about the phys- 
ical content of a possible zero in t-matrix in the near- 
threshold regime translates one-to-one also to production 
reactions. However, this is not correct for there is, in ad- 
dition to the hadronic production, also the production via 
the quark component possible. This piece can be expressed 
via the spectral density. Indeed, as seen from the right 
panel Fig. [51 the corresponding production amplitude is 
given by: 

Mq^~FqGqo{E)tqh{k,E), ^^E, (52) 

where is the production amplitude for the bare quark 
state by the point-like source, Gqo{E) = l/{Eo~E) is 
the quark state free Green's function, and tqh(k,E) is the 
t-matrix element responsible for the quark-meson transi- 
tion. From the corresponding Lippmann-Schwinger equa- 
tion. 



tqh{k,E)= f{k) 



f{p)t{p,k,E) 



(53) 



and with the help of Eqs. ([T 
found in the form: 



13, the latter can be 



tqhik,E) = J d?pxk{p)f{p) ^{E- Eo)ck{E). (54) 

Thus, the near-threshold production rate via the quark 
component of the wave function is simply 



dBrq{MiM2) 
dE 



— const X w{E), 



(55) 



where u;(i?) is defined in Eq. ([20)1 . Then, using definition 
(PO)) and the explicit form of the coefficient Ck{E), which 
follows from Eq. p^ . 



Ck{E) 



9f ^ 1 
Stt^P E-Ec 



(56) 



w{E) 



1 kgf {Ef-Ecf 
27r\VF\^ {E-EcY 



(57) 



Obviously, this implies that dBrq/dE does not employ a 
zero. 

In reality one expects both mentioned mechanisms to 
contribute to the production of the physical resonance, 
with a relative importance depending on the production 
mechanism. Therefore, we may write for the full produc- 
tion rate: 



dBr(MiM2) 
dE 



const X k 



\E-Ec + r{Ef -Ec) 
\Vf\^{E-EcY 



(58) 

where r is a real number containing, among other contri- 
butions, the ratio of the production rates via the quark 
state and the hadronic state. A priori no estimate of r 
is possible thus, even if there were a system that shows a 
zero in the scattering amplitude near threshold and allow- 
ing for all the conclusions of the previous section, this zero 
might well be shielded in the production reaction. On the 
other hand, li Ec ^ A it requires a delicate fine tuning 
of the parameter r to produce a zero in the amplitude of 
Eq. (l58)) . Thus, if there is a zero observed in a production 
amplitude, it at least suggests the presence of both quark 
states as well as a hadronic molecule. 

Note that the zero in the lineshape of the X(3872) 
predicted to occur in a B-decay amplitude reported in 
Ref. [231124] is of a different kind, for it originates as a 
coupled-channel effect of two nearby continuum channels. 
In this work, on the other hand, the effects in the presence 
of only a single continuum channel were discussed. 



6 Implications for lattice QCD 

From the previous discussions it should be clear that the 
presence of a zero in the scattering amplitude contains im- 
portant information about the system studied. However, 
this information is no longer visible in the production am- 
plitude, which is the quantity accessible experimentally. 

Fortunately scattering observables, in particular phase 
shifts, are accessible in lattice QCD from a study of volume 
dependencies of the energy spectra |25| . In simulations of 
full QCD both quark states as well as molecular states are 
present and will influence all correlators, as long as some 
overlap exists with the interpolating fields used. Thus, the 
results of those simulations need to be interpreted in the 
same way as experimental results call for an interpreta- 
tion, if one wants to gain some insight into the dynami- 
cal mechanisms that lead to the structure formation. One 
method pushed recently is to work with a large basis of 
interpolating fields and to use the resulting eigenvectors of 
the correlator matrices to interpret the findings [251[?7] . In 
Ref. [22] the scattering length is shown to be an important 
quantity also to extract information about the nature of 
states from lattice studies. The discussion above provides 
one with an additional approach: in order to understand 
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the interplay of quark and meson degrees of freedom one 
may investigate the near-threshold phase shifts. If a zero 
in the i-matrix, corresponding to a zero in the phase shifts, 
occurs, it immediately signals the presence of both quark 
poles as well as potential scattering poles with a relatively 
weak coupling as described in Sec. HI 

To illustrate this point we discuss briefly the phase 
shift (its cotangent) for Case (i) introduced above, 



k cot S = — 



(59) 



This expression employs a pole at the position of the t- 
matrix zero, as shown in Fig. [S] We expect that such a 
structure could be extracted from studies within lattice 
QCD. 



7 Conclusions 

In this paper we have presented a dynamical scheme in 
which a near-threshold zero in the mesonic t-matrix ap- 
pears as a consequence of the interplay between quark 
states and a nonperturbatively interacting hadron-hadron 
continuum. The appearance of such a zero invalidates the 
effective-range expansion and corresponds to three near- 
threshold poles of the t-matrix. The effect on hadronic 
observables in elastic scattering is proved to be drastic. 

However, a near-threshold t-matrix zero could exist 
only if several requirements are met. First, one needs the 
direct interaction in the mesonic channel to be strong 
enough to support a bound or virtual state. Second, a 
nearby bare quark state should exist, with a weak cou- 
pling to the mesonic channel. While, in principle, such a 
situation cannot be excluded, it is highly accidental. 

Without such special arrangements the effective-range 
formulae are valid, and the conclusions of Refs. pTllT^fT^ 
[HlfTSlfTB] hold true: a large and negative effective range 
corresponds to a compact quark state, while a small ef- 
fective range means that the state is composite. In the 
former case there are two near-threshold poles in the t- 
matrix, while in the latter case there is only one pole. 

In the two-pole case one can definitely state that the 
resonance is generated by an s-channel diagram. The cou- 
pling of this quark state to the mesonic continuum is not 
large. In the one-pole situation no model-independent in- 
sight into the underlying dynamics is possible and one can 
only indicate that the resonance is generated dynamically 
but, in the language of meson exchange, the binding po- 
tential could emerge from either s-channel of i-channel 
exchanges — or a mixture of both. It seems not pos- 
sible to decide between these scenarios model indepen- 
dently. The three pole scenario, on the other hand, calls 
for the presence of both s-channel interactions as well 
as potential scattering — traditionally identified with t- 
channel exchanges. This situation can be identified easily 
by the break-down of the effective-range expansion. Al- 
though scattering experiments are not possible for most 
unstable particles, an alternative access to scattering ob- 
servables could be provided by studies of lattice QCD. 



We also discussed the implications for production reac- 
tions. Unfortunately it turns out that for the lineshapes of 
production reactions no statement is possible model inde- 
pendently based solely on information from the scattering 
amplitudes. 

The work was supported in parts by funds provided from the 
Helmholtz Association (grants VH-NG-222, VH-VI-231), by 
the DFG (grants SFB/TR 16 and 436 RUS 113/991/0-1), by 
the EU HadronPhysics2 project, by the RFFI (grants RFFI-09- 
02-91342-NNIOa and RFFI-09-02-00629a), and by the Presi- 
dential programme for support of the leading scientific schools 
(grants NSh-4961.2008.2 and NSh-4568.2008.2). The work of 
V.B., Yu.K., A.K., and A.N. was supported by the State Cor- 
poration of Russian Federation "Rosatom". A. N. would like 
to acknowledge the support of the non-profit "Dynasty" foun- 
dation and ICFPM and of the grant PTDC/FIS/70843/2006- 
Fisica. 



References 

1. E. S. Swanson, Phys. Rept. 429 (2006) 243. 

2. E. Klempt and A. Zaitsev, Phys. Rept. 454 (2007) 1. 

3. M. B. Voloshin, Prog. Part. Nucl. Phys. 61 (2008) 455. 

4. S. Godfrey and S. L. Olsen, Ann. Rev. Nucl. Part. Sci. 58 
(2008) 51. 

5. E. Braaten, arXiv:0808.2948 [hep-ph]. 

6. S.-K. Choi et al. [Belle Collaboration], Phys. Rev. Lett. 91 
(2003) 262001. 

7. B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 
95 (2005) 142001. 

8. B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 
98 (2007) 212001. 

9. S.-K. Choi et al. [Belle Collaboration], Phys. Rev. Lett. 
100 (2008) 142001. 

10. R. Mizuk et al. [Belle CoUaboration], Phys. Rev. D 78 
(2008) 072004. 

11. S. Weinberg, Phys. Rev. 130 (1963) 776. 

12. S. Weinberg, Phys. Rev. 131 (1963) 440. 

13. S. Weinberg, Phys. Rev. 137 (1965) B672. 

14. V. Baru et al, Phys. Lett. B 586 (2004) 53. 

15. V. Baru et al, Eur. Phys. Journ. A 23 (2005) 523. 

16. D. Morgan, Nucl. Phys. A 543 (1992) 632. 

17. L. N. Bogdanova, G. M. Hale, and V. E. Markushin, Phys. 
Rev. C 44 (1991) 1289. 

18. M. Doring, C. Hanhart, F. Huang, S. Krewald and U. G. 
Meifiner, Phys. Lett. B 681 (2009) 26. 

19. M. Doring, C. Hanhart, F. Huang, S. Krewald and U. G. 
Meifiner, Nucl. Phys. A 829 (2009) 170. 

20. S. Flatte, Phys. Lett. B 63 (1976) 224. 

21. C. Hanhart, J. R. Pelaez, and G. Rios, Phys. Rev. Lett. 
100 (2008) 152001. 

22. F. K. Guo, C. Hanhart, and U. G. Meifiner, Eur. Phys. J. 
A 40 (2009) 171. 

23. E. Braaten and M. Lu, Phys. Rev. D 77 (2008) 014029. 

24. M. B. Voloshin, Phys. Rev. D 76 (2007) 014007. 

25. M. Liischer, Commun. Math. Phys. 105 (1986) 153. 

26. J. J. Dudek, R. G. Edwards, M. J. Peardon, D. G. 
Richards, and C. E. Thomas, Phys. Rev. Lett. 103 (2009) 
262001. 

27. J. J. Dudek, R. G. Edwards, M. Mathur, and D. G. 
Richards, Phys. Rev. D 77 (2008) 034501. 



V. Baru et ai: Interplay of quark and meson degrees of freedom in a near-threshold resonance 



5000 r 



4000 r 



3000 - 



2000 



<j [mb] 





• E [MeV] 



E [MeV] 



Fig. 1. Case (i). Upper panel: the pole structure in the fc-plane (left plot) and the Z-factor — see Eq. (|36p (right plot) versus 
the coupling constant g/. Lower panel: the elastic scattering cross section (left plot) and the spectral density (right plot) versus 
the energy for gj — 0.01 (solid line), gr/ = 0.03 (dashed line), and gj = 0.1 (dotted line). 




Fig. 2. The poles motion in the s-plane for Case (i) at the first Riemann sheet (left plot) and on the second Riemann sheet 
(right plot). The unitarity cut is depicted in grey. 
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Fig. 3. Case (ii). Upper panel: the pole structure in the fc-plane. Lower panel: the elastic scattering cross section (left plot) and 
the spectral density (right plot) versus the energy for <?/ = 0.01 (solid line), g/ = 0.03 (dashed line), and = 0.1 (dotted line). 
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Fig. 4. Case (iii-a). Upper panel: the pole structure in the fc-plane (left plot) and the Z- factors — see Eq. p6|l (right plot) 
versus the coupling constant gf. Lower panel: the elastic scattering cross section (left plot) and the spectral density (right plot) 
versus the energy for gf = 0.001 (solid line), gf = 0.02 (dashed line), and (?/ = 0.3 (dotted line). 
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Fig. 5. Case (iii-b). Upper panel: the pole structure in the fc-plane (left plot) and the Z-factors (right plot) versus the coupling 
constant gf. Lower panel: the elastic scattering cross section (left plot) and the spectral density (right plot) versus the energy 
for gf — 0.001 (solid line), gf = 0.02 (dashed line), and gf = 0.3 (dotted line). 
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Fig. 6. Case (iv-a). Upper panel: the pole structure in the fc-plane (left plot) and the Z-factor — see Eq. p6|) (right plot) 
versus the coupling constant gf. Lower panel: the elastic scattering cross section (left plot) and the spectral density (right plot) 
versus the energy for gf = 0.01 (solid line), gf = 0.13 (dashed line), gf = 0.3 (dotted line), and g = 0.5 (dash-dotted line). 
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Fig. 7. Case (iv-b). Upper panel: the pole structure in the fc-plane (left plot) and the Z-factor — see Eq. (|36|l (right plot) 
versus the coupling constant gf. Lower panel: the elastic scattering cross section (left plot) and the spectral density (right plot) 
versus the energy for gf = 0.01 (solid line), gj — 0.07 (dashed line), gf — 0.15 (dotted line), and g — 0.5 (dash-dotted line). 
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Fig. 8. Diagrams for the two meson production via the hadronic component (left plot) and via the quark component (right 
plot). 




Fig. 9. The phase shift for Case (i) for gf — 0.01 (solid line), gf — 0.02 (dashed line), and gf = 0.1 (dotted line). 



